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There are no reasons why the energy spectra of the relic gravitons amplified by the pumping 
action of the background geometry should not increase at high frequencies. A typical example of this 
behavior are quintessential inflationary models where the slopes of the energy spectra can be either 
blue or mildly violet. In comparing the predictions of scenarios leading to blue and violet graviton 
spectra we face the problem of correctly deriving the sensitivities of the interferometric detectors. 
Indeed the expression of the signal-to-noise ratio not only depends upon the noise power spectra 
of the detectors but also upon the spectral form of the signal and, therefore, one can reasonably 
expect that models with different spectral behaviors will produce different signal-to-noise ratios. 
^\ ' By assuming monotonic (blue) spectra of relic gravitons we will give general expressions for the 

0^ , signal-to-noise ratio in this class of models. As an example we studied the case of quintessential 

■ gravitons. The minimum achievable sensitivity to h\ JIgw of different pairs of detectors is computed, 

and compared with the theoretical expectations. 
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I. INTRODUCTION 



Gravitational wave astronomy, experimental cosmology and high energy physics will soon experience a boost thanks 
to the forthcoming interferometric detectors. From a theoretical point of view it is then interesting to compare our 
theoretical expectations/speculations with the foreseen sensitivities of the various devices in a frequency range which 
■ complements and greatly extends the information we can derive from the analysis of the microwave sky and of its 
CN| temperature fluctuations. 

By focusing our attention on relic gravitons of primordial origin we can say that virtually every variation in the time 
evolution of the curvature scale can imprint important informations on the stochastic gravitational wave background 
0. The problem is that the precise evolution of the curvature scale is not known. Different cosmological scenarios, 
based on different physical models of the early Universe may lead to different energy spectra of relic gravitons and 
£^ ] this crucial theoretical indetermination can affect the expected signal. 

Of particular interest seems to be the case where the logarithmic energy density of the relic gravitons (in critical 
units) grows [^,0in the frequency region explored by the interferometric detectors (i.e., approximately between few 
X • Hz and 10 kHz) In this range we can parametrize the energy density of the relic gravitons pew a t the present 



time 770 as 



n GW (f, Vo ) = -^j = n(%)?(/,»jd) (1.1) 



where VL(rjo) denotes the typical amplitude of the spectrum and q(f, rjo) is a monotonic function of the frequency at 
least in the interval 1 Hz < / < 10 kHz. Both D,(r)o) and q(f, rjo) can depend on the parameters of the particular 
model. The assumption that q(f,r]o) is monotonic can certainly be seen as a restriction of our analysis, but, at the 
same time, we can notice that the models with growing logarithmic energy spectra which wer e discussed up to now 



in the litterature fit in our choice for q(f, t]q). Within the parametrization defined in Eq. (1.1) we will be discussing 
the cases where the spectral slope a (i.e., a — dq(f,i] )/df) is either blue (i.e., < a < 1) or violet (i.e., a > 1). In 
general we could have also the case a < (red spectra) and a — (flat spectrum) . Flat spectra have been extensively 
studied in the context of ordinary inflationary models j^] and in relation to cosmic string models ]l0| . 

Blue and violet spectra are physically peculiar since they are typically produced in models which are different from 
the ones leading to flat spectra. In quintessential inflationary models [ p"l| the logarithmic energy spectra are typically 
blue jl2| . This is due to the fact that in this class of models an ordinary inflationary phase is followed by an expanding 
phase whose dynamics is driven by an effective equation of state which is stiffer than radiation Jl3[ ] . Since the equation 
of state (after the end of inflation) is stiffer than the one of radiation, then the Universe will expand slower than in a 
radiation dominated phase and, therefore, a turns out to be at most one (up to logarithmic corrections). 
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In string cosmological models jl4| the graviton spectra can be either blue (if the physical scale corresponding to a 
present frequency of 100 Hz went out of the horizon during the string phase) or violet (if the relevant scale crossed 
the horizon during the dilaton driven phase). 

The purpose of this paper is to analyze the sensitivity of pairs of interferometric detectors to blue and violet 
spectra of relic quintessential gravitons. The reason for such an exercise is twofold. On one hand violet and blue 
spectra, owing to their growth in frequency, might provide signals which are larger than in the case of flat inflationary 
spectra. On the other hand the sensitivity to blue spectra from quintessential inflation can be different from the one 
computed in the case of flat spectra from ordinary inflationary models. Indeed, it is sometimes common practice to 
compare the theoretical energy density of the produced gravitons with the sensitivity of various interferometers to a 
flat spectrum. This is, strictly speaking, arbitrary even if, sometimes this procedure might lead to correct order of 
magnitude estimates. 

In order to illustrate qualitatively this point let us consider the general expression of the signal-to-noise ratio (SNR) 
in the case of correlation of two detectors of arbitrary geometry for an observation time T. By assuming that the 
intrinsic noises of the detectors are stationary, gaussian, uncorrelated, much larger in amplitude than the gravitational 
strain, and statistically independent on the strain itself, one has JTH] [l^] :[] 

SNR * = J0- 2 fvt{ f d/ jf«5swm r, (1 . 2) 

2V2ir 2 (Jo /»si 1) (/)ai 2 >(/)/ 

(Ho is the present value of the Hubble parameter and F depends upon the geometry of the two detectors; in the 
case of the correlation between two interferometers F = 2/5). In Eq. (1.2), S„(f) is the (one-sided) noise power 



spectrum of the fc-th (fc = 1,2) detector, while 7(/) is the overlap reduction function [|17|,[18| which is determined by 
the relative locations and orientations of the two detectors. This function cuts off (effectively) the integrand at a 
frequency / ~ l/ 2d, where d is the separation between the two detectors. 



From eq. (1.2) we can see that the frequency dependence of the signal directly enters in the determination of the 
SNR and, therefore, we can expect different values of the integral depending upon the relative frequency dependence 
of the signal and of the noise power spectra associated with the detectors. Hence, in order to get precise information 
on the sensitivities of various detectors to blue and violet spectra we have to evaluate the SNR for each specific model 
at hand. 

The analysis of the SNR is certainly compelling if we want to confront quantitatively our theoretical conclusions with 
the forthcoming data. Owing to the differrence among the various logarithmic energy spectra of the relic gravitons we 
can wonder if different detector pairs can be more or less sensitive to a specific theoretical model. We will try, when 
possible, to state our conclusions in such a way that our results could be used not only in the specific cases discussed 
in the present paper. 

The plan of the paper is the following. In Section II we will review the basic features of blue spectra arising in 
quintessential inflationary models. In Section III we will set up the basic definitions and conventions concerning the 
evaluation of the SNR. In Section IV we will be mainly concerned with the analysis of the achievable sensitivities to 
some specific theoretical model. Section V contains our concluding remarks. 

II. BLUE AND VIOLET GRAVITON SPECTRA 
A. Basic bounds 

Blue and violet logarithmic energy spectra of relic gravitons are phenomenologically allowed |19| . At low frequencies 



the most constraining bound |2C|] comes from the COBE observations 21 1 of the first (thirty) multipole moments of the 
temperature fluctuations in the microwave sky which implies that hg ^gw(/o ; Vo) has to be smaller than 6.9 x 10 -11 



1 Notice that, with this definition, the SNR turns out to be the square root of the one used in refs. JL5|-|l8|. The reason for our 
definition lies in the remark that the cross-correlation between the outputs si,2(£) of the detectors is defined as: 



S = / ' dt [ 7 dt' Sl (t)s 2 (t')Q(t,t'), 

J-T/2 J -T/2 



rT/2 f T/2 

dt 

-T/2 J-T/2 

where Q is a filter function. Since S is quadratic in the signals, with the usual definitions, it contributes to the SNR squared. 
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for frequencies of the order of Hq. At intermediate frequencies (i.e., f p ~ 10 8 Hz) the pulsar timing measurements 
p2f imply that riGw(/pj77fi) should not exceed 10 -8 . In order to be compatible with the homogeneous and isotropic 



nucleosynthesis scenario J23 24 we should require that 



K J GW (/,?7o)dln/ < 0.2n 7 (r? )^ =s 5 x 10" 6 , (2.1) 

where f2 7 (7/o) = 2.6 x 1 0~ 5 h^ 2 is the fraction of critical energy density in the form of radiation at the present 
observation time. In Eq. (2.1) the integral extends over all the modes present inside the horizon at the nucleosynthesis 
time. In the case of blue and violet logarithmic energy spectra the COBE and pulsar bounds are less relevant than the 
nucleosynthesis one and it is certainly allowed to have growing spectra without conflicting with any of the bounds.^] 

B. Quintessential Spectra 

Recent measurements of the red-shift luminosity relation in type la supernovae p7| suggest the presence of an 
effective cosmological term whose energy density can be as large as 0.8 in critical units. Needless to say that this energy 
density is huge if compared with cosmological constant one would guess, for instance, from electroweak (spontaneous) 
symmetry breaking, i.e., p\ ~ (250 GeV) 4 . In order to cope with this problem various models have been proposed |2q | 
and some of them rely on the existence of some scalar field (the quintessence field) whose effective potential has no 
minimum |p9| . Therefore, according to this proposal the evolution of the quintessence field is dominated today by the 
potential providing then the wanted (time-dependent) cosmological term. In the past the evolution of the quintessence 
field is in general not dominated by the potential. The crucial idea behind quintessential inflationary models is the 
identification of the inflaton <fi with the quintessence field Therefore, the inflaton/quintessence potential V{4>) 

will lead to a slow-rolling phase of de Sitter type for <f> < and it will have no minimum for <fi > 0. Hence, after 
the inflationary epoch (but prior to nucleosynthesis ) the Universe will be dominated by <j) 2 . This means, physically, 
that the effective speed of sound of the sources driving the background geometry during the post-inflationary phase 
will be drastically different from the one of radiation (i.e., c s = l/\/3 in natual units) and it will have a typical stiff 
form (i.e., c s — 1). The fact that in the post-inflationary phase the effective speed of sound equals the speed of light 
has important implications for the gravitational wave spectra as it was investigated in the past for a broad range of 
equations of state stiffer than radiation (i.e., l/\/3 < c s < 1) |l^]. The conclusion is that if an inflationary phase is 
followed by a phase whose effective equation of state is stiffer than radiation, then, the high frequency branch of the 
graviton spectra will grow in frequency. The tilt depends upon the speed of sound and it is, in our notations, 

We can immediately see that for all the range of stiff equations of state (i.e., l/\/3 < c s < 1) we will have that 
< a < 1 . The case a = corresponds to c s = l/\/3- This simply means that if the inflationary phase is 
immediately followed by the ordinary radiation-dominated phase the spectrum will be (as we know very well) flat. 
The case c s = 1 is the most interesting for the case of quintessential inflation. In this case the tilt is maximal (i.e., 
a = 1). Moreover, a more precise calculation jl2|,[l3| shows that the graviton spectrum is indeed logarithmically 
corrected as 

?(/, r*>) = £ln 2 (£). (2.3) 

It is amusing to notice that this logarithmic correction occurs only in the case c s — 1 bu t not in the case of the other 
stiff (post-inflationary) background. The typical frequency /i(?7o) appearing in Eq. J2.3| ) is given, today, by 

/i(r? ) = 1132 A/7 3 ( ^ ) GHz. (2.4) 

gth 



2 Notice that the nucleosynthesis bound refers to the case where the underlying nucleosynthesis model is ho mogen eous and 
isotropic. The presence of magnetic fields and/or matter-antimatter fluctuations can slightly alter the picture |25[|2q] . 
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Apart from the dependence upon the number of relativistic degrees of freedom (i.e., gdcc = 3.36 and <?th = 106.75 ) 
which is a trivial consequence of the red-sfhit, /i(t?o) does also depend upon N s which is the number of (minimally 
coupled) scalar degrees of freedom present during the inflationary phase. The amplitude of the spectrum depends 
upon N s as 

— 1 64 x 1CT 5 

nfo) = Ng — ■ (2.5) 

The reason for the presence of N s is that all the minimally coupled scalar degrees of freedom present during the 
inflationary phase will be amplified sharing approximately the same spectrum of the two polarizations of the gravitons 
pO| . The main physical difference is that the N s scalars are directly coupled to fermions and, therefore, they will decay 
and thermalize thanks to gauge interactions p3] |. If minimally coupled scalars would not be present (i.e., N s = 0) the 
model w ould not be consistent since the Universe will be dominated by gravitons with (non-thermal) spectrum given 
by Eq. ( [2.3] ). The energy density of the quanta associated with the minimally coupled scalars, amplified thanks to 
the background transition from the inflationary phase to the stiff phase, will decrease with the Universe expansion 
as a~ A whereas the energy density of the background will decrease as a -6 . The moment at which the energy density 
of the background becomes sub-leading marks the beginning of the radiation dominated phase and it takes place at 
a (present) frequency of the order of the mH z |13| . Notice that this frequency has been obtained by requiring the 
reheating mechanism to be only gravitational ]30|. This assumption migh be relaxed by considering different reheating 
mechanisms [[nj (see also |}2| ) . In order to satisfy the nucleosynthesis constraint in the framework of a quintessential 
model with gravitational reheating |30| we have to demand that pl| , p^ ] 

o / \ 1/3 

4-[—) <°- 07 > ( 2 - 6 ) 

where the factor of 3 counts the two polarizations of the gravitons but also the quanta associated with the inflaton 
and g n = 10.75 is the number of spin degrees of freedom at t n . For frequencies /(r/ ) > /i(r?o) the spectra of the 
produced gravitons are exponentially suppressed as exp [— f/fi]. This is a general feature of the spectra of massless 
particles produced thanks to the pumping action of the background geometry jl2|][] 



III. SIGNAL-TO-NOISE RATIO FOR MONOTONIC BLUE SPECTRA 



In order to detect a stochastic gravitational wave background in an optimal way we have to correlate the outputs of 
two (or more) detectors |l5|-|l8| . The signal received by a single detector can be thought as the sum of two components: 
the signal (given by the stochastic background itself) and the noise associated with each detector's measurement. The 
noise level associated with a single detector is, in general, larger than the expected theoretical signal. This statement 
holds for most of the single (operating and/or foreseen) gravitational waves detectors (with the possible exception of 
the LISA space interferometer ||). Suppose now that instead of a single detectors we have a couple of detectors or, 
ideally, a network of detectors. The signal registered at each detector will be 

s i = h i (t) + n i (t), (3.1) 

where the index i labels each different detector. If the detectors are sufficiently far apart the ensamble average of the 
Fourier components of the noises is stochastically distributed which means that 

«(/K(/')> = ^(/-.n^ ) (i/i), (3.2) 

where SVi(|/|) is the one-sided noise power spectrum which is usually expressed in seconds. The very same quantity 
can be defined for the signal. By then assuming the noise levels to be statistically independent of the gravitational 
strain registered by the detectors we obtain Eq. ( |l.2f ). 

Consider now the case of two correlated interferometers and define the following rescalcd quantities: 



3 Quintessential graviton spectra have, in general, three branches: a soft branch (for 1CP 18 Hz <i / <, 1CP 16 Hz), a semi-hard 
branch (for 10 -16 Hz <, f <! 1CP 3 Hz) and a hard branch which is the one mainly discussed in the present paper. The reason 
for this choice is obvious since the noise power spectra of the interferometric detectors are defined in a band which falls in the 
region of the hard branch of the theoretical spectrum. 
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. E« = S 7 { l) /So (i = 1,2); 

• v = ///o; 

• «gw(/) - n(/o)w(/). 

(in this Section we will not write the explicit dependence of the theoretical quantities upon rj : they are meant to be 
considered at the present time). Notice that / is (approximately) the frequency where the noise power spectra are 
minimal and r2cw(/o) is the graviton (logarithmic) energy density at the frequency /o- Therefore the signal-to-noise 
ratio can be expressed as: 

SNR 2 = J*L VT J , (3.3) 

where we defined the (dimension-less) integral 

J2 r dy ±m:im . ( 3. 4) 

A ^"(/M^W) 
From this last expression we can deduce that the minimum detectable h\ r2cw(/o) is given by (1 yr = tt x 10 7 s) 

/?5/2 q / \ 1/2 

,2 o ( t.\ A n ^ If 32 I yr > omo2 

For example, by taking /o = 100 Hz and So = 10 -44 s, we get 



h 2 fi GW (/o) - 4.0 x 10 32 I ) SNR (3.5) 



^n G w(100Hz) c 4 -° X J 1 °' 7 (^) V2 SNR 2 . (3.6) 

Therefore, the estimate of the sensitivity of cross-correlation measurements between two detectors to a given spectrum 
Hgw(/) reduces, in our case, to the calculation of the integral J defined in Eq. (3.4). Given a specific theoretical 
spectrum, J can be numerically determined for the wanted pair of detectors. 



IV. ACHIEVABLE SENSITIVITIES FOR QUINTESSENTIAL SPECTRA 



Consider first the case of the two LIGO detectors (located at Hanford, WA and Livingston, LA) in their "advanced" 
versions. From the knowledge of the geographical locations and orientations of these detectors |53|, the overlap 
reduction function can be calculated |l^|l8|], and the result is reported in Fig. 2 of ref. ]l8fl . As function of the 
frequency, 7 has its first ze ro a t 64 Hz and it falls rapidly at higher frequency. This behavior allows to restrict the 
integration domain in Eq. (ETj) to the region / < 10 kHz (i.e., v < 100). We assumed identical construction of the 



(i) 



two detectors (i.e., Sn 
ref. |3J, namely (see Fig. 



H) 



) . For the rescaled noise power spectrum of each detector we used the analytical fit of 



£„(/) = 



3G 



fa 



hi ^ 4, 



I (r/a 



a j2 

• \3 J 



f<fb 

h<f<Y 

f > F/a 



(4.1) 



with 



1.96 x 10~ 2 T = 1.6 f a = 68 Hz f b = 10 Hz. 
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FIG. 1. We report the rescaled noise power spectra of the LIGO and VIRGO detectors used for the calculation of the 
signal-to-noise ratio. 

In the case of a flat spectrum (i.e., a — 0, tu(f) = 1) we find J ~ 6.1 x 10 3 , which implies 

1 :\ 1 / 2 



/i5ft G w(100Hz) ~ 6.5 x 10" 11 y^fr J SNR ( 4 - 2 ) 

in close agreement with the estimate obtained in ref. jTsj. 

The minimum detectable flaw f° r quintessential gravitons can be obtained by recalling that 



<"(/) = rs- * - 

In i^i \ v\ 



For /o = 100 Hz, numerical integration gives: 



1 Q3 , ,„ 

J ~ —= — {6.91 + 21.36 lni/i + 26.52 In 2 i/ x + 15.68 In 3 i/ x + 3.78 In 4 v x } 
In v\ 



J - foon u^r^2 ftW) (4-3) 



or, taking into account Eq. fl2.4|) , in terms of N s : 

1.6 x 10 7 
3.0-lni\g 

with 



P^(iV s ) ~ 1.07 - 4.62 x 10~ 2 lniV s + 7.52 x 10~ 4 In 2 iV s 
-5.44 x 10~ 6 ln 3 AT s + 1.48 x 10~ 8 ln 4 7V s . 



By inserting this expression in Eq. ( p.6| ), one has: 



^(WdOOBh) c 2.5 x ID- 4 (88 ^" ( y )!> (^) 1/2 SNR 2 (4.4) 



By assuming for N s the minimum value compatible with Eq. ( |2.6[ ) (i.e., AT S = 21), wc obtain: 

/ign GW (100Hz) ~ 1.8 x 10~ 10 ( J SNR 2 (4.5) 
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As we can see by comparing Eq. ( |4.2D with Eq. (4.5) the minimum detectable f2cw(100 Hz) is slightly larger 
for growing spectra. This a general result that is simply related to the structure of J. For the special value of N s 
considered this difference is roughly of a factor of 2. Another important point to stress is that for both the graviton 
spectra considered, as a consequence of the frequency behavior of 7(/) and the presence of the weighing factor v~ e 
in the integrand, the main contribution to the integral J comes from the region / < 100 Hz. The cut-off introduced 
by the overlap reduction function is not so relevant: by assuming j(f) — 1 over the whole integration domain (i.e., 
considering the correlation of one of the detector with itself), the sensitivity increases only by a factor 2.4 in the case 
of a flat spectrum, and 3.6 in the case of the quintessential one. This means that the only way to get a substantial 
rise in sensitivity lies in the improvement of the noise characteristics of the detectors in the low-frequency region. 

As a comparison we considered also the sensitivity that could be obtained at VIRGO in the (purely hypothetical) 
case in which the detector now under construction at Cascina, near Pisa (Italy), were correlated with a second 
interferometer located at about 50 km from the first and with the same orientation.^] The overlap reduction function 
for this correlation has its first zero at a frequency / ~ 3 kHz (see Fig. 0). 




f (Hz) 



FIG. 2. The overlap reduction function for the correlation of VIRGO with a coaligned interferometer whose central (corner) 
station is located at (43.2 N, 10.9 E), d ~ 58.0 km (Italy). 



Also in this case we assumed that the detectors are identical and for the common rescaled noise power spectrum 
we used the analytical parametrization given in ref. p5[ (see Fig. |l|) 



where 



£n(/) 



s 2 J -4 



1+4- 



f a = 500 Hz 



ft = 2Hz 



/ < fb 
f>fb 



(4.6) 



51 = 3.46 x 10~ 6 

5 2 = 6.60 x 10~ 2 
£ 3 = 3.24 x 10~ 2 



4 For illustrative purposes, we assumed, within our example, 50 km as the minimum distance sufficient to decorrelate local 
seismic and e.m. noises. This hypothesis might be proven to be correct and it is certainly justified in the spirit of this exercise. 
However, at the moment, we do not have any indication either against or in favor of our assumption. 
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In the c ase of flat spectrum, limiting the numerical integration to 10 kHz, we obtain J ~ 5.5 and, therefore, according 
to Eq. (U) 



ft,o«Gw(100 Hz) ~ 7.2 x 10" 

In the case of quintessential inflation, for fo = 100 Hz, we have: 
1 



1 yr 

T 



1/2 



SNR 



(4.7) 



J 



\n 2 v x 



{ 5.79 - 0.30 In v x + 31.20 In 2 v x + 6.11 In 3 v x + 12.91 In 4 v x } 



1/2 



or, in terms of N s 



with 



1.6 x 10 4 
■8.0-lniV s ) : 



Pv(N s ) 



(4.8) 



P$(N S ) ~ 3.10 - 0.14 ln7V s + 2.37 x 10~ 3 In 2 N s 



17.84 x 10" 6 ln 3 iV s + 5.04 x 10~ 8 In 4 N s 



Therefore, from Eq. ([03) one has: 



/i 2 ^Gw(100Hz) ~ 2.5 x 10 



_ u (88.0-lniV s ) 2 / lyr\ 1/2 



Pv(N s ) 



T 



SNR 



that, for N s = 21 gives: 



/M<> ( .rA (HK)llz) - 1.1 « 10" 7 ( ^ 



1/2 



SNR 



At a frequency of 100 Hz the theoretical signal can be expressed as 



/i^n GW (100Hz) = 7V S " 3/4 x 10" 15 



2220.07 - 50.46 In N s + 0.281n 2 iV s 



(4.9) 



(4.10) 



(4.11) 



as a function of N s . In Fig. || this function (full thin line) is compared with the s ens itivity of LIGO-WA*LIGO-LA 
and VIRGO*VIRGO (full thick lines) obtained from, respectively, Eqs. |u|) and (gj), assuming T = 1 yr and SNR 
= 1. We can clearly see that our signal is always below the achievable sensitivities. Notice that, if we assume purely 
gravitational reheating N s > 21. 

One could think that, thanks to the sharp growth of the spectrum, the signal could be strong enough around 10 
kHz, namely at the extreme border of the interferometers band. Indeed around fo — 10 kHz, the theoretical signal is 
given by: 



hgfiGwCLOkHz) = A7 3/4 x 10~ 15 



1387.81 - 39.89 In N 8 + 0.28 In 2 N s 



(4.12) 



We see that the situation does not change qualitatively. In fact it is certainly true that around 10 kHz the signal is 
larger but the sensitivity is also smaller. In fact, repeating the calculation for fo — 10 kHz, in the case N s = 21 we 
obtain: 



hi ^Gw(lOkHz) 



1.1 x 10~ 8 
6.7 x 10~ 6 



1 yr 

1 yr 

T 



1/2 



SNR 



1/2 



SNR 



LIGO - WA * LIGO - LA 
VIRGO * VIRGO 



(4.13) 



If we compare Eqs. (4.13) with Eqs. d4.5| ) and ( 4. 10] ) we see that the minimum detectable signal gets larger the 
larger is the spectral frequency. Therefore, the mismatch appearent from Fig. ^| between the theoretical signal and 
the experimental sensitivity will remain practically unchanged. 
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FIG. 3. We report the theoretical amplitude computed in Eq. (4.11) (full thin line) and the associated sensitivities computed 
in Eqs. (O) and (O) for T = 1 yr and SNR = 1 (full thick lines). 



V. CONCLUDING REMARKS 



In this paper we precisely computed the sensitivity of pairs of interferometric detectors to blue and mildly violet 
spectra of relic gravitons. Our investigation can be of general relevance for any model predicting non flat spectra 
of relic gravitons. We analyzed the correlation of the two LIGO detectors in their "advanced" phase. On a more 
speculative ground we investigated the theoretical possibility of the correlation of VIRGO with an identical, coaligned, 
interferometer located very near to it. 

As a test for our techniques we first discussed the case of a flat spectrum which has been discussed in the past. 
We then applied our results to the case of quintessential inflationary models whose graviton spectra are, in general, 
characterized by three "branches". A soft branch (in the far infra-red of the graviton spectrum around 10 -18 -10 -16 
Hz), a semi-hard branch (between 10~ 16 Hz and 10 -3 Hz) and a truly hard branch ranging, approximately, from 10~ 3 
Hz to 100 GHz. Since the interferometers band is located, roughly, between few Hz and 10 kHz, the relevant signal 
will come from the hard branch of the spectrum whose associated energy density appears in the signal-to-noise ratio 
with blue (or mildly violet) slope. In the hard branch the energy density of quintessential gravitons is maximal for 
frequencies in the range of the GHz. In this region Hq f^cw can be as large as 10 -6 . In spite of the fact quintessential 
spectra are growing in frequency the predicted signal is still too small and below the sensitivity achievable by the 
advanced LIGO detectors. The reason for the smallness of the signal in the region / ~ 1 kHz is twofold. On one 
hand we have to enforce the nucleosynthesis bound on the spectrum. On the other hand, because of the gravitational 
reheating mechanism adopted, the number of (minimally coupled) scalar degrees of freedom needs to be large. It 
might be possible, in principle, that different reheating mechanisms could change the signal for frequencies comparable 
with the window of the interferometers. Therefore, the analysis presented in this paper seems to suggest that new 
techniques (possibly based on electromagnetic detectors p6[) operating in the GHz region should be used in order to 
directly detect quintessential gravitons. 
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